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A  study  of  the  contact  mechanics  of  a  probe  tip  interacting  with  a  quartz  crystal  microbalance 
(QCM)  has  been  performed,  involving  simultaneous  measurements  of  normal  load,  displacement, 
and  contact  stiffness  with  changes  in  QCM  resonant  frequency.  For  metal-metal  and  glass-metal 
contacts  in  air,  the  QCM  frequency  shifts  were  observed  to  be  positive,  and  directly  proportional  to 
the  contact  area  as  inferred  from  the  contact  stiffness.  Interfacial  characteristics  of  the  probe-tip 
contact  (elasticity,  contact  size,  and  an  estimate  of  the  number  of  contacting  asperities)  were 
deduced  by  extending  a  prior  model  of  single  asperity  contact  to  the  case  of  multiple  contacts.  The 
extended  model  clarifies  a  number  of  seemingly  disparate  experimental  results  that  have  been 
reported  in  the  literature.  ©  2001  American  Institute  of  Physics.  [DOI:  10.1063/1.1413493] 


I.  INTRODUCTION 

Knowledge  of  the  physical  properties  of  small  contact¬ 
ing  asperities  is  a  topic  that  is  of  increasing  importance  to  a 
number  of  macroscopic  and  microscopic  applications.  Nu¬ 
merous  experimental  characterizations  of  the  physical  prop¬ 
erties  of  nanoscale  contacts  have  been  performed  by  means 
of  nanoindentation, 1-4  atomic  force  microscopy5-8  and  scan¬ 
ning  tunneling  microscopy  (STM).9'10  More  recently,  it  has 
become  recognized  that  a  probe  tip-quartz  crystal  microbal¬ 
ance  (QCM)11-20  or  closely  related  geometries21  could  also 
have  enormous  potential  as  a  sensitive  probe  of  interfacial 
physical  properties.  The  principle  underlying  such  experi¬ 
ments  is  that  the  extremely  sharp  resonance  of  a  quartz  crys¬ 
tal  is  sensitive  to  interaction  forces  applied  to  its  surfaces. 
Changes  in  resonance  frequency  /0  and  quality  factor  Q  re¬ 
flect  the  additional  energy  storage  and  loss  accompanying  the 
interaction.  In  general,  all  four  combinations  of  increasing  or 
decreasing  /0  and  increasing  or  decreasing  Q  are  possible.22 
Beginning  in  1994,  the  frequency  shift  of  a  QCM  in  response 
to  a  contacting  probe  tip  was  reported  by  groups  indepen¬ 
dently  seeking  to  characterize  asperity  contacts.  Krim  and 
co-workers1112  reported  negative  frequency  shifts  and  de¬ 
creases  in  Q  due  to  the  presence  of  a  STM  tip  for  adsorbed 
water  films  of  varying  thickness.  Sasaki  et  al. 13-14  reported 
the  use  of  STM  for  “scanning  shear  stress  microscopy,”  at¬ 
tributing  the  observed  QCM  frequency  shifts  to  be  directly 
proportional  to  interfacial  shear  stress.  The  spatial  contrast  in 
these  and  similar  experiments1516  has  been  variously  attrib¬ 
uted  to  local  friction,  elasticity,  or  viscoelasticity. 

Two  groups  have  recently  evaluated  the  QCM-probe  in¬ 
teraction  to  interpret  the  physical  properties  of  such  contacts. 
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Laschitsch  and  Johannsmann 17  19  (LJ)  studied  a  number  of 
materials  combinations  and  reported  positive  frequency 
shifts  and  decreased  quality  factors.  They  presented  a  model, 
assuming  a  single  asperity  contact,  which  predicted  positive 
frequency  shifts  linearly  proportional  to  contact  radius.  In 
contrast,  Flanigan  et  al.20  directly  measured  the  contact  area 
between  a  QCM  and  a  polymer  gel  cap  and  observed  a  nega¬ 
tive  frequency  shift  directly  proportional  to  the  contact  area. 
A  model  based  on  known  expressions  for  the  QCM  fre¬ 
quency  shift  on  full  immersion  in  a  fluid  medium  was  used  to 
deduce  the  viscoelastic  properties  of  the  gels. 

The  seemingly  disparate  results  reported  in  the  literature 
stem  largely  from  the  fact  that  no  experimental  group  has 
performed  measurements  in  completely  characterized  situa¬ 
tions  (normal  force,  contact  area,  etc.),  and  thus  untested 
assumptions  have  necessarily  been  incorporated  into  the 
analyses.  In  order  to  simultaneously  measure  the  relevant 
physical  parameters  to  interpret  such  measurements,  we  have 
used  a  combined  nanoindenter-QCM  apparatus.  This  allows 
independent  measurement  of  the  applied  load,  normal  dis¬ 
placement,  and  contact  stiffness,  while  simultaneously  moni¬ 
toring  the  frequency  shift  of  the  QCM.  The  advantage  of 
measuring  the  displacement  and  contact  stiffness  throughout 
the  QCM-probe  interaction  is  that  it  allows  quantitative 
evaluation  of  the  contact  mechanics  governing  the  interac¬ 
tion.  For  the  probe  and  surface  configuration  employed  in 
our  experiments,  we  observed  that  the  frequency  shift  of  the 
QCM  was  positive  and  linearly  proportional  to  the  contact 
area.  Interfacial  characteristics  of  the  probe-tip  contact 
(elasticity,  contact  size,  and  an  estimate  of  the  number  of 
contacting  asperities)  were  deduced  by  extending  the  LJ 
model  of  single  asperity  contact  to  the  case  of  multiple  con¬ 
tacts.  The  extended  model  clarifies  a  number  of  the  experi¬ 
mental  results  that  have  been  reported  in  the  literature. 
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FIG.  1.  Schematic  diagram  of  the  experimental  apparatus. 


II.  EXPERIMENT 

Figure  1  shows  a  diagram  of  the  experimental  setup 
combining  a  depth-sensing  nanoindenter  and  a  QCM.  The 
nanoindenter  consists  of  a  load-displacement  transducer  with 
electrostatic  force  actuation  and  displacement  sensing  elec¬ 
tronics  (Hysitron,  Inc.,  Minneapolis,  MN).  With  a  custom 
force  modulation  setup,23  this  instrument  enables  dynamic 
measurements  (at  about  105  Hz  in  these  experiments)  of  the 
contact  stiffness  normal  to  the  substrate,  Kn ,  and  avoids 
problems  associated  with  thermal  drift.  It  is  well  established 
that  for  a  spherical  indenter  tip,  Kn  is  directly  proportional  to 
the  contact  radius,  ac,  and  given  by 

Kn  =  2E*ac ,  (1) 

where  E*  is  the  reduced  Young’s  modulus.24  With  known 
reduced  modulus  of  the  materials  in  contact,  the  contact  ra¬ 
dius  or  area  can  be  directly  calculated  from  the  measured 
contact  stiffness.2 

The  QCM  was  clamped  in  a  holder  that  secures  the  crys¬ 
tal  at  three  points  along  its  edge  and  mounted  on  an  XYZ 
piezoscanner  base  (Digital  Instruments,  Santa  Barbara,  CA). 
The  crystals  used  were  fundamental  mode  5  MHz  AT  (trans¬ 
verse  shear  mode)-cut  circular  quartz  blanks  with  1  cm  di¬ 
ameter.  Silver  electrodes  ~100  nm  thick  and  0.6  cm  in  di¬ 
ameter  were  thermally  evaporated  onto  the  blanks.  The 
crystals  were  driven  with  a  Clapp  feedback  oscillator 
circuit,25  and  the  resonant  frequency  and  amplitude26  were 
monitored  with  a  frequency  counter  and  an  oscilloscope,  re¬ 
spectively.  The  quality  factors  of  the  resonators  were  mea¬ 
sured  with  a  simple  “ring-down”  technique27  and  were  near 
10.  Two  probes  were  used,  a  fused  silica  sphere  (Duke  Sci¬ 
entific,  Palo  Alto,  CA),  0.5  mm  radius,  and  an  etched  tung¬ 
sten  wire,  ~12  /xm  tip  radius.  The  experiments  were  con¬ 
ducted  at  atmospheric  pressure  in  a  dry  nitrogen  environment 
with  1%-10%  relative  humidity. 

During  a  typical  experiment,  the  QCM  was  first  brought 
into  contact  with  the  indenter  tip  (~1  /xN  load  setpoint)  us¬ 
ing  the  piezoscanner.  Care  was  taken  to  center  the  indenter 
tip  on  the  QCM  electrode  to  maximize  QCM  sensitivity. 
Then,  a  preprogrammed  indentation  cycle  (load/unload)  was 
initiated.  The  normal  force  was  ramped  in  either  a  staircase 
fashion  with  plateaus  lasting  30  s,  or  with  a  linear  ramp  at 
loading  rates  of  25  to  50  /xN/s  with  a  10-15  s  hold  before 
unloading;  maximum  loads  were  100  to  500  /xN.  The  data 
acquisition  program  simultaneously  recorded  the  normal 
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FIG.  2.  Plot  of  QCM  frequency  shift  vs  normal  load  for  contacts  with  the 
silica  sphere  (0.5  mm  radius  of  curvature)  and  W  tip  (radius  —12  pm). 


load,  displacement  of  the  probe,  contact  stiffness,  and  QCM 
resonant  frequency.  For  the  staircase  ramp  experiments,  the 
resolution  of  the  frequency  shift  was  ±0.02  ppm  (±0.1  Hz), 
using  a  gate  time  of  1  s.  For  the  linear  ramps,  a  faster  gate 
time  of  0.1  s  was  used,  which  reduced  the  frequency  shift 
resolution  to  0.5  Hz.  The  attractive  forces  during  the  initial 
contact  of  these  probe  tips  with  the  QCM  surface  were  in¬ 
sufficient  to  cause  observable  frequency  shifts,  as  measured 
using  a  force-distance  curve  approach  method.28  In  these 
experiments,  a  minimum  load  of  approximately  5  /xN  was 
needed  to  observe  a  +0.1  Hz  frequency  shift  with  the  12  pm 
W  tip. 


III.  EXPERIMENTAL  RESULTS 


Figure  2  shows  an  example  of  normal  load  versus  QCM 
frequency  shift  for  both  the  W  tip  and  silica  sphere  contacts. 
The  frequency  response  of  the  QCM  was  linearly  propor¬ 
tional  to  applied  load  over  a  range  up  to  500  /xN.  For  these 
contacts,  loads  of  30-50  /xN  were  needed  to  obtain  a  1  Hz 
frequency  shift  on  the  QCM.  Figure  3  shows  normal  load 
versus  contact  stiffness  data  for  the  silica  sphere.  The  re¬ 
sponse  is  clearly  nonlinear,  and  is  best  fit  with  a  square-root 
dependency.  Deviation  at  the  higher  loads  is  the  result  of 
error  in  the  measurement  from  a  decrease  in  displacement 
amplitude  with  increasing  contact  stiffness.  Similar  behavior 
was  observed  for  the  contacts  made  using  the  W  tip  (not 
shown),  which  were  also  best  fit  with  a  square-root  depen¬ 
dency. 

The  contact  area  can  be  estimated  from  the  contact  stiff¬ 
ness  measurements  using  Eq.  (1).  The  reduced  modulus  is 
given  by 


(2) 


where  £,  and  v,  are  the  Young’s  modulus  and  Poisson’s 
ratios  for  the  substrate  and  the  probe.  In  calculating  the  re¬ 
duced  modulus,  we  have  taken  the  thin  Ag  film  on  quartz  as 
a  single  system  with  modulus  70  GPa  and  Poisson’s  ratio 
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FIG.  3.  Contact  stiffness  vs  normal  load  for  the  silica  sphere.  The  response 
is  best  fit  with  a  square-root  dependency. 


FIG.  5.  Frequency  shift  vs  contact  area  for  the  silica  sphere  and  W  tip 
during  loading  to  500  /jN. 


0.2. 29  We  assumed  a  modulus  of  70  GPa  and  345  GPa  and 
Poisson’s  ratio  of  0.2  and  0.28,  for  the  silica  sphere  and  W 
tip,  respectively. 30  The  resulting  reduced  moduli  were  36 
GPa  for  the  silica  sphere  and  61  GPa  for  the  W  tip. 

Figure  4  shows  a  plot  of  the  calculated  contact  area  ver¬ 
sus  normal  load  for  both  the  silica  sphere  and  W  tip.  Linear 
fits  are  shown,  with  R 2  values  of  0.99.  The  linear  relationship 
between  contact  area  and  normal  load  demonstrates  that  the 
contact  is  non-Hertzian.  [For  a  sphere  against  flat  geometry 
in  Flertzian  (elastic)  contact,  the  contact  area  varies  with  load 
as  L2/3.]  The  linear  relationship  is  however  consistent  with  an 
interface  that  is  multicontact  and/or  undergoing  plastic  defor¬ 
mation.  Either  of  these  conditions  alone  is  sufficient  to  give 
the  observed  proportionality;31  and  both  are  expected  to  play 
a  role.  Load-displacement  data  (not  shown)  from  indenta¬ 
tions  on  the  silvered  region  of  the  QCM  showed  evidence  of 
plastic  deformation  under  these  loading  conditions.  Further¬ 
more,  the  contact  areas  reported  here  always  span  a  range 
much  larger  than  the  observed  asperities  on  the  QCM  elec¬ 
trode  surface,32  so  multicontact  interfaces  are  expected. 


FIG.  4.  Contact  area  vs  normal  load  calculated  from  the  contact  stiffness 
data  during  loading  for  the  silica  sphere  and  W  tip. 


Figure  5  shows  a  plot  comparing  QCM  frequency  shift 
versus  the  contact  area  for  both  the  W  and  silica  contacts 
during  loading  to  500  /jN.  For  both  contacts,  the  relationship 
is  linear  ( R 2  values  0.98-0.99),  with  a  slope  unique  to  each 
interface.  (We  note  that  departures  from  the  linear  relation¬ 
ship  of  Fig.  5  were  sometimes  observed  for  higher  loads, 
indicating  that  the  contact  may  be  entering  a  new  regime  of 
deformation  and/or  wear.) 

The  linear  relationship  between  frequency  shift  and  con¬ 
tact  area  shown  in  Fig.  5  is  an  apparent  contradiction  with 
the  LJ17  model,  which  predicts  a  linear  proportionality  be¬ 
tween  frequency  shift  and  contact  radius.  On  the  other  hand, 
the  linear  relationship  between  contact  area  and  frequency 
shift  is  consistent  with  the  results  of  Flanigan  et  al.,20  but  the 
frequency  shift  is  in  the  opposite  direction.  We  propose  that 
the  multicontact  nature  of  many  probe -QCM  interfaces  is  an 
important  consideration.  We  now  show  that  our  results  can 
be  accounted  for  by  modifying  the  LJ  model  accordingly.  In 
doing  so,  the  proportionality  constants  apparent  in  Fig.  5 
may  be  analyzed  in  terms  of  the  number  and  size  of  contacts 
at  the  interface. 


IV.  MODELING 

In  an  earlier  article,  LJ  considered  a  geometry  in  which 
the  radius  of  the  contact  zone  was  much  smaller  than  the 
wavelength  of  shear  waves  propagating  into  the  probe.17  In 
this  model,  the  frequency  shift  A /  and  “loss  tangent” 
AT/ A/  are  each  linearly  proportional  to  the  contact  radius 
and  can  be  identified  with  an  “added  stiffness”  term.  (Here, 
T  is  the  half  bandwidth  at  half  maximum  power).  However, 
the  derived  proportionality  between  resonance  shifts  and 
contact  radius  could  not  be  verified,  as  this  required  an  ex¬ 
plicit  measurement  of  the  contact  size  and  its  dependence  on 
experimental  parameters.  The  shear  stress  acting  at  the 
probe-QCM  interface  can  be  characterized  using  a  wave 
model  for  the  propagation  of  acoustic  waves  from  the  QCM 
into  the  probe.  LJ  employed  this  approach  to  derive  the 
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FIG.  6.  The  BVD  equivalent  circuit  for  a  QCM,  with  a  small  perturbing 
load  impedance  ZL.  L  is  proportional  to  the  mass  per  unit  area  of  the  quartz 
blank,  1/C  is  proportional  to  the  stiffness  per  unit  area  G/d  (where  G  is  the 
shear  modulus  and  d  the  thickness  of  the  quartz),  and  R  is  proportional  to  the 
viscous  dissipation  inside  the  quartz  blank. 


changes  in  QCM  resonance  parameters  using  a  well-known 
equivalent  circuit  technique,  which  we  briefly  describe  in  the 
following  paragraph. 

Within  linear  circuit  analysis,  a  QCM  can  be  described 
as  a  damped  harmonic  oscillator  by  the  lumped-element 
Butterworth-van  Dyke  (BVD)  equivalent  circuit,  shown  in 
Fig.  6.  The  BVD  circuit  is  derived  from  the  more  rigorous 
transmission  line  model  due  to  Mason  with  several  reason¬ 
able  approximations.33  To  be  treated  in  this  context,  acoustic 
interactions  involving  the  probe  and  quartz  crystal  are  as¬ 
sumed  to  be  both  linear  in  displacement  and  one¬ 
dimensional,  i.e.,  the  acoustic  waves  are  planar  with  no  lat¬ 
eral  variations.34  The  latter  assumption  is  certainly  not  true  in 
practice  since  the  shear  mode  amplitude  of  a  QCM  decays  in 
a  Gaussian  manner  from  the  center  of  the  electrode 
outwards,35  and  the  interaction  with  the  probe  tip  is  very 
localized.  (Normalization  to  the  electrode  active  area  intro¬ 
duced  later  on  allows  this  otherwise  one-dimensional  model 
to  account  for  much  of  the  lateral  variation.)  Nevertheless, 
we  state  these  assumptions  and  indicate  that  their  limits  of 
applicability  must  be  tested. 

The  acoustic  interaction  between  the  probe  and  QCM 
was  treated  by  LJ  as  an  added  load  impedance  ZL  in  series 
with  the  other  elements  of  the  acoustic  branch.  For  semi¬ 
infinite  media,  ZL  is  simply  the  mechanical  impedance 
evaluated  at  the  interacting  surface  of  the  QCM.  The  me¬ 
chanical  impedance  is  defined  by 

ZL  =  vlv,  (3) 


where  cr  is  the  applied  shear  stress  and  v  is  the  surface  ve¬ 
locity.  ZL  is  a  complex  quantity  and  must  be  normalized  to 
the  electrode  area.  Since  the  changes  in  resonance  param¬ 
eters  are  only  a  few  parts  per  million,  the  interaction  between 
the  probe  and  QCM  can  be  treated  as  a  small  perturbation. 
Within  the  BVD  damped  harmonic  oscillator  model,  the 
changes  in  resonant  frequency  and  half  bandwidth  are  given, 
to  first  order,  by: 


A  /= 


and 


/o  Im[Zj 
77  Zq 


(4) 


/o  R e[Z£] 

77  Zq 


(5) 


where  f0  is  the  fundamental  frequency  and  Zq  =  8.8 
X 106  kg  m  2  s  1  is  the  acoustic  impedance  of  AT-cut 


FIG.  7.  Schematic  diagram  of  the  probe-QCM  interface. 


quartz.  Notice  that  the  loss  tangent  AT/ A/  is  simply  the 
negative  ratio  of  the  real  and  imaginary  parts  of  ZL .  For  light 
damping,  the  change  in  dissipation  D  is  given  by 


AD  =  A 


1  \  I  AE\  l2AT\ 
Q]=A\2^Ej=\y^)' 


(6) 


where  A E  is  the  energy  dissipated  per  cycle  and  E  is  the 
energy  stored  in  the  system. 

To  evaluate  the  mechanical  impedance  and  determine  the 
change  in  resonance  parameters,  one  must  propose  an  acous¬ 
tic  model  and  calculate  the  shear  stress  at  the  interface.  We 
can  show  that  a  “near-field”  acoustic  geometry  applies  to 
our  system,  similar  to  that  described  by  LJ.  The  wavelength 
of  5  MHz  acoustic  waves  in  all  materials  used  here  is  near 
600  /rm,36  whereas  the  contact  radii  are  between  50  and  300 
nm.  Therefore,  the  contact  radius  is  at  least  three  orders  of 
magnitude  smaller  than  the  wavelength  of  sound.  As  such, 
we  expect  sound  to  radiate  into  the  tip  as  spherical  waves 
decaying  as  et r,  where  e  is  a  characteristic  size  of  the  true 
contact  regions  (asperities)  through  which  sound  propagates, 
and  r  is  the  radial  distance  from  the  contact.  The  simplest 
way  to  represent  the  acoustic  deformation  of  the  contact  is 


u(r,t)  =  U  o 


e 

r+  e 


ei(wt-kr ) 


(7) 


where  u0  is  the  amplitude  of  the  QCM  shear  oscillation,  &>  is 
the  angular  frequency,  and  k  is  the  wave  number.  As  shown 
in  Fig.  7,  e  equals  the  contact  radius  ac  for  a  single-contact 
interface  where  the  number  of  contacts  N=  1 . 

For  a  multicontact  interface,  we  propose  e  is  approxi¬ 
mately  constant  as  both  N  and  the  contact  area  increase  in 
proportion  to  the  load.37  In  this  case,  e  is  determined  by  the 
nature  of  the  roughness  on  the  scale  of  the  contact  zone. 
Evaluating  the  shear  stress  at  the  interface,  r  =  0,  we  find: 


du 
K  — 
dr 


Ku, 


r= 0 


(8) 


where  K  is  an  effective  shear  modulus  of  the  interface.  The 
surface  mechanical  impedance  can  then  be  evaluated,  taking 
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care  to  normalize  to  the  area  of  shear  waves  in  the  QCM  by 

2  2 

including  a  factor  ac/ae,  where  ae  is  the  radius  of  the  active 
area  (essentially  the  electrode  area). 


a 

Zr= -4 


a2  du l dt 


r  =  Q 


a]  K 

■  —  — ( ke-i ). 
we 


Using  Eqs.  (4)  and  (5)  we  find 


A  /= 


fo  Ka2c 


irZq(oae 


a  r= 


/  o 


“7 - 2  Kka'c 

n TZqa)ae 


(9) 

(10) 

(ID 


and 


(12) 


The  interface  is  characterized  by  K,  a  c ,  and  e.  From  a 
dimensional  analysis,  the  quantity  Ka2/  e  defines  a  spring 
constant  for  the  increased  stiffness  of  the  system  and  positive 
frequency  shift.  The  increase  in  half  bandwidth  and  dissipa¬ 
tion  due  to  sound  propagation  depends  on  the  total  contact 
area  independent  of  the  roughness.  More  importantly,  we 
find  that  for  interfaces  where  e-constant,  the  frequency  shift 
is  proportional  to  the  contact  area  while  the  loss  tangent  re¬ 
mains  constant.  On  the  other  hand,  in  the  limit  of  single¬ 
asperity  interfaces  where  e=ac,  we  recover  the  result  of  LJ, 
i.e.,  both  the  frequency  shift  and  loss  tangent  are  propor¬ 
tional  to  contact  radius. 

While  LJ  had  no  independent  means  of  measuring  the 
contact  size,  the  data  published  in  Ref.  17  allow  one  to  de¬ 
termine  whether  the  dependence  of  the  frequency  shift  and 
loss  tangent  on  the  contact  radius  is  self-consistent.  For  the 
smooth  lubricated  contact,  A /  and  AT/ A/  are  observed  to 
increase  with  contact  radius  in  the  same  fashion,  agreeing 
with  their  model.  In  contrast,  for  rougher  metal-metal  con¬ 
tacts,  A /  increases  by  at  least  a  factor  of  4  while  AI7A / 
remains  essentially  constant.  This  apparent  departure  from 
their  model  is  explained  within  the  roughness  picture  pre¬ 
sented  here.38 

We  may  analyze  the  slopes  of  the  data  presented  in  Fig. 
5  to  show  that  our  model  gives  reasonable  results  for  the  size 
and  number  of  contacting  regions.  According  to  Eq.  (10),  for 
oscillations  in  the  fundamental  mode  u>  =  2  irf0,  the  slope  is 
given  by 

A  /  K 

A~c=  2^-ZqAee’  (13) 

where  Ae=7ra2,  Ae  =  2.8XlO-5  m2,  and  Zq=  8.8 
X  106  kg  m  2  s  '.  It  is  an  approximation  to  consider  only 
one  value  for  the  elasticity  of  the  interface  in  our  model, 
especially  given  the  proposed  spherical  geometry  of  the 
wave  fronts  in  the  probe.  Since  the  QCM  response  is  corre¬ 
lated  with  the  lateral  deflection  of  the  contact,  one  choice  for 
K  is  the  reduced  shear  modulus,  G*,  of  the  interface:39 


1  2— iq  2  —  v2 

G*~  Gx  +  G2 


(14) 


TABLE  I.  Calculation  of  contact  region  (asperity)  radius  for  the  W  and 
silica  probes. 


Probe  material 

A f/Ac  (from  Fig.  5) 

Asperity  contact  radius,  e 

W  tip  (R=  12  fim) 

102±1  Hz/jttm2 

26  nm 

Silica  sphere  (R= 0.5  mm) 

16.6±0.3  Hz//zm2 

99  nm 

Again,  we  consider  the  Ag  film  and  quartz  to  be  a  single 
system,  with  shear  modulus  of  29  GPa,  while  the  W  tip  and 
silica  sphere  have  shear  moduli  of  135  GPa  and  29  GPa, 
respectively.  For  the  corresponding  interfaces,  the  reduced 
shear  moduli  are  then  13  and  8  GPa.  Fitting  Eq.  (13)  to  the 
measured  slopes  in  Fig.  5  allows  us  to  estimate  of  the  aver¬ 
age  size  of  asperity  contacts,  shown  in  Table  I.  These  values 
are  reasonable  since  asperities  in  this  size  range  are  com¬ 
monly  observed  in  STM  images  of  QCM  electrode  surfaces 
(see,  for  example.  Refs.  27  and  32).  The  range  of  areas  dis¬ 
played  in  Fig.  5  corresponds  to  a  range  of  up  to  330  nm  in 
contact  radius.  Our  data  and  model  are  both  consistent  with 
an  interface  involving  a  collection  of  distinct  contact  regions 
with  a  constant  average  size  e.  Considering  e  to  be  the  radius 
of  asperities,  the  number  of  asperity  contacts  may  be  esti¬ 
mated  from  the  area  ratio:  n  =  a2cl e2.  For  the  W  tip,  we  find 
the  estimated  number  of  asperity  contacts  varies  from  1  to 
~70  as  contact  area  increases  in  Fig.  5.  For  the  silica  sphere, 
the  number  varies  from  1  to  —10.  Again,  these  values  seem 
reasonable  for  our  simple  model.  Presumably,  one  could  de¬ 
termine  the  overall  size  of  the  contact  zone  by  estimating  the 
separation  between  contacts.  However,  this  would  require  a 
more  detailed  analysis  of  the  asperity  height  distribution  and 
degree  of  deformation  than  the  model  justifies. 

The  model  presented  here  indicates  that  changes  in  fre¬ 
quency  shift  are  simply  proportional  to  changes  in  local  elas¬ 
ticity  as  long  as  the  contact  area  and  roughness  remain  un¬ 
changed.  This  clarifies  the  mechanism  of  lateral  contrast  in 
local  elasticity  imaging  setups  based  on  scanning  probe- 
QCM  combinations. 13-16  This  model  does  not  contradict  the 
observation  of  negative  frequency  shifts  in  experiments  with 
low  modulus  polymer  probes  on  QCM  surfaces.20  One  can 
show  that  the  lower  shear  moduli  and  significantly  larger 
contact  areas  imply  that  the  wavelength  of  sound  in  the  poly¬ 
mer  is  not  larger  than  the  contact  radius,  so  the  interaction  is 
no  longer  in  the  near-field  acoustic  limit.  In  this  case,  theo¬ 
ries  based  on  plane  waves  apply,  which  predict  the  observed 
negative  frequency  shifts.22  This  might  have  been  the  case  as 
well  for  the  negative  frequency  shifts  reported  in  Ref.  1 1  for 
tip-QCM  interactions  in  the  presence  of  water  films. 

The  model  we  have  described  assumes  the  contact  re¬ 
gions  do  not  slip.  Equation  (8)  shows  that  for  a  shear  modu¬ 
lus  of  about  10  GPa,  QCM  oscillation  amplitude  u0~  10  nm 
(a  typical  value  for  the  QCM’s  in  these  experiments),27  and 
e— 100  nm,  the  shear  stresses  applied  to  the  interface  may 
extend  into  the  GPa  range,  enough  to  cause  significant  plas¬ 
tic  deformation  and  wear.  Indeed,  a  transfer  film  was  ob¬ 
served,  apparently  of  silver  from  the  electrode  surface,  on  the 
silica  sphere  after  the  experiments.  Indentation  curves  ob¬ 
tained  with  the  QCM  oscillating  and  static  were  compared, 
and  it  was  found  that  the  depth  of  penetration  of  the  tip  into 
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the  oscillating  QCM  was  nearly  twice  that  of  the  static  QCM, 
although  both  exhibited  hysteresis.  This  suggests  that  in  ad¬ 
dition  to  elastoplastic  deformation,  wear  was  indeed  occur¬ 
ring  at  some  point  during  the  indentations,  implying  fretting 
or  slip.  It  is  well  known  that  the  onset  of  slip  depends  on 
both  oscillation  amplitude  and  the  load,  with  slip  favored  at 
larger  amplitudes  and  smaller  loads.24  However,  the  results 
presented  here  are  not  sufficient  to  determine  the  degree  of 
slip  at  the  interfaces,  and  the  model  describes  our  data  quite 
well  in  terms  of  no-slip  boundary  conditions. 

Nevertheless,  both  mass-and-spring  mechanical  models 
and  a  recently  published  physical  model22  indicate  that  the 
onset  of  interfacial  slip  is  associated  with  a  decrease  in  the 
magnitude  of  the  frequency  shift  and  a  broadening  of  the 
resonance  due  to  frictional  energy  losses.  These  consider¬ 
ations  indicate  that  it  should  be  possible  to  investigate  inter¬ 
facial  shear  strengths,  fretting,  and  high-speed  sliding  fric¬ 
tion  using  probe-QCM  techniques.  The  pressures  and 
sliding  speeds  accessed  would  be  relevant  to  the  operation  of 
microelectromechanical  systems  as  well  as  macroscopic  de¬ 
vices,  for  the  former  of  which  there  historically  have  been 
very  few  nano-  or  microscale  techniques  operating  in  the 
appropriate  physical  regime.40 

Finally,  we  see  from  our  results  that  it  should  be  possible 
to  carry  out  nanomechanical  studies  with  probe-QCM  setups 
in  which  the  number  of  contact  regions  is  controllably  varied 
over  the  mesoscopic  range  from  a  single  to  hundreds  of  con¬ 
tacts.  Such  studies  would  benefit  greatly  from  specially  pre¬ 
pared  probe  and/or  substrate  surfaces  (with  desirable  mor¬ 
phologies  or  material  properties),  increased  frequency  shift 
sensitivity,  and  a  wide  variation  of  the  QCM  vibrational  am¬ 
plitude,  for  the  possibility  of  studying  both  nonsliding  and 
sliding  contacts. 

V.  CONCLUSIONS 

We  have  investigated  the  contact  mechanics  occurring  at 
the  interface  between  a  probe  tip  and  QCM  surface.  For  the 
tip-surface  combinations  chosen,  the  data  and  modeling  pre¬ 
sented  are  consistent  with  an  interface  involving  many  dis¬ 
crete  contact  regions.  The  change  in  resonant  frequency  of 
the  QCM  is  directly  proportional  to  the  contact  area  and 
reflects  the  elasticity  of  the  interface  and  the  number  and 
average  size  of  contact  regions. 
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